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1. INTRODUCTION 

Environmental charging bf geosynchronous spacecraft to potehtials of 
thousands Of volts has been experimentally observed. Previous attempts to model 
Spacecraft charging have employed techniques which are limited to simplified 
geometries artd Symmetry assumptions. In this paper, wc describe a new compu- 
tational model which Can simulate the charging Of complex geometrical objects in 
three dimensions. We present two sample calculations. In the first problem, the 
capacitance to infinity of a complex object similar to a Satellite with Solar array 
paddles is calculated, the Second problem concerns the dynamical charging of a 
conducting cube partially covered with a thin dielectric film, tn this calculation, 
the phOtoemlSsiOn results in differential charging of the object. 


2. THEORV 

the interaction of a satellite artd the magnetosphere can be separated into two 
parts, the first is the particle deposition, charge transport, and electrical 
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properties associated with charged particles Iniplngihg upon the satellite, "The 
eecohd part Is the sell- consistent artiblent and photbplasma Interactions with the 
electric field. The field must Satisfy boundary conditions on the satellite consist- 
ent with the charge deposited oh the satellite, to solve both parts of this problem .. 
completely and self- consistently for general, ambient plasmas IS a formidable 
task. Here, we shall be concerned only with a limited (albeit very important) 
range of plasma environments (a hot magnetosphere). As a reSult, certain appro- 
xlmations reduce thfe magnitude of the problem. 

The timescales of phenomena which Occur on a Spacecraft in the magneto- 
sphere range from nanoseconds to hours (Table 1). The lower end of this range 
is associated with effects such as the discharging of electrical circuit elements 
and electromagnetic waVe phenomena. At the upper end Of the range, slow vari- 
ations in the magnetospheric environment are important. In this. paper, we shall 
be concerned with the intermediate timescale range, from milliseconds to seconds. 
This range is determined by the charging time of the surface of a spacecraft by 
magnetospheric electron currents. 


Table 1. Characteristic Times for Charged Spacecraft in the Magnetosphere 


Phenomenon 

Time 

Breakdown in circuit elements 

Tg ~ lo"® - lo“® sec 

Charging of bare conducting surface 

T 10 sec 

c 

Differential charging of thin dielectric 
overlylttg conductor 

~ 1 sec 

Charge rediStributiofi In a dielectric 


Change in environmental cdndltibne 

^ 1 to ^10^ sec 


Before entering upon the analysis of these phenomena, it is useful to set the 
scale of the various prbceSSes involved. These are listed in Table 2, and, in each 
case, the treatments which must be applied to describe the relevant field and 
particle phenomena are indicated, in the magnetosphere, the plasma Sheath sur- 
rounding a spacecraft requires a particle description. This Is necessary ort 
account of the very long mean free paths and long Debye lengths which occur 
in these hot, diftuse plasmas, felectromagnetlc treatments are needed only for 
describing effects such as transient surface current phenomena resUitlng from 
arcing excitations. Particle dynamics must be followed If sheath plasma 
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Table 2. CoUisioblees SatelUte-Plasnfia Sheath Mbdele, L ta the Spacecraft 
Stae, top ie the-PlaSrtia Fr'equency and ta the Debye Length 


Model 

Fields 

Particle 

Treatrtient 

Timescale 

1. Electromagnetic 

Maxwell's equations 

Dynamic 

L'c 

— lO”® sec 

2. Quasistatic 

Poisson's equations 
(L « CWp* 

Dynamic 

^ 10 ^ sec 

3. Equilibrium 

{a) Poisson 

(b) Laplace 
(Xjj » L) 

Static 

>10*^ sec 


oscillatory behavior is important. ^ On a longer timescale, the plasma is 
characterized by an equilibrium-particle distribution. ^ This is the range which 
is considered here. Finally, at the longest timescale, the bf havior is determined 
by changes in the environment Or redistribution of charges within dielectrics. 

We Shall nO^ describe the considerations underlying our analysis of the inter- 
mediate timescale phenomena. Let uS consider a spacecraft with a Spherical 
conducting Surface. With a radius R (cm) and a charging current density 
j (A/cm^). the time taken to charge the spac.e.fi,r5ft..lQ.A. potential V is 

C V 


where is the capacitance of the sjiacecraft With respect to infinity, and ie 
given by 

C = R esu 
00 

i 10*^^ R F . 

With the follovving values 

R = 100 crn 
V = 1 q3 V 

J = 0. 5 X lO’® A/cm^ 
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the chei^glng time Is 




■i 2 X 10"® see 


However, an most spacecraft, isrge ax'eaS Sre hot here conductors, but are 
covered by thin, insulating dielectrics overlylrtg conducting substrates. lA these 
cases, the capacitance of the dielectric. Cp. Is important. leather than the capaci- 
tance with respect to infinity. The dielectric capacitance is 


jid esu/cm2 


10 


-12 


4>rd 


F/cm^ 


With k thickness of 40 mil fd ^0. 1 cm), the charging time Tj^ is new 

-0 


10 V 

’’n '■ 

D 4ird j 

1. 6 Set 


10 


L2 XO. 5 X 10 




The voltage buildup between a Conductor and a dielectric insulator thus occurs 
very much more slowly than the buildup on a bare conductor. Asa result, differ- 
ential charging of a spacecraft takes place on a timescale longer than charging. 
Later, we shall describe the development of such a phenomenon o^fet the time- 
scale range < t < 

Undet* the conditions found In magnetospheric substorms, essential simplifies' 
tions can be made in the modeling of the ck "ging. in particular. We demonstrate 
below for hot. low density plasmas where 

» L. 

is the Debye length and L a characteristic object dimension, that if surface 
fkotentials on the satellite are of the order of the filasma temperature, one makes 
only a very small error by neglecting the ambient Space charge density in 
Poisson's etiuation. lliis approximation, when justified, greatly reduces the 
amount of computation necessary to determine satellite boteAtials. 
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I^et us ekattitne tlt$ effect of a lOrge ambient charge density fluctuation ih a 
d ° 10 keV, ® 10 cm*^ magrtetos^herio plasma, The Debye length of such a 
plasma is given by 

‘‘n ^ "^43 / — i 2 X 10^ cm 

“ \l “e 

1 = 200 m. 

ft has been shown^ that in equilibrium plasmas, maximum charge variations are 
of Order of the ambient charge density. It follows then that the magnitude of a 
potential associated With a spherical charge density fluctuation Of 1 m in radius 
is at most 

♦ ~ »r^ n^e esu 

= |#‘l0^ X 10 X 4. 8 X 10’ 

~ 2 X 10“^ StatVoltS 
~ R X 10"^ Volts 

which is several orders of magnitude less than the satellite surface potentials. 

Another useful quantity to examine iS the relative amount Of charge on a 
sphere of a meter radius charged up to the ambient temperature to the amount 
Of Space charge such a volume would contain. The surface charge on a sphere of 
radius r is 

**surface ‘ ■ 


The Space charge in such a plaSrna is 


plasma “ 3 "e* 


The ratio of these two charges is 


**plasma 
**Surface ^ 


= » » r n 


J=— _ * 

T 


< 10 


•5 


323 


I 

I 

I 


ThuH, wt* have strong reasohs to believe that the gross potential features 
RvuToundln'! ad object whose dirtlenslons are much smaller than a Debye length 
oikI whose surface potentials are comparable to the plasma temperature can be 
ealculntcd Ignoring ambient (as distinct from photosheath) space charge effects 


:l. M MKllU.U. KS 

The dynamical model consists of twc parts, namely: (1) the calculation of 
surface charge densities and net, d’.arging currents, given a potential distribution, 
and (2) the calculation of the potential subject to free space and appropriate satel- 
lite boundaiy conditions. Brief descriptions of the techniques used are given 
below. Further details will be presented at a later date. 


;l. I Surface Churcinp (iuleulutiun 

We require the incident and outgoing currents respectively at sur- 

face points ~r * r^. The net charging current is then 

^net ' ^in^o^ ’ ^out^*’o^ 
where 

1, (r ) “ /* d^ V ^ f. (v , r ) 

••in' o J o o in o' o 

i Ar ) - f V V f . (v , r ) 

■'out o J o o out o o 

The distribution f^^^ is assumed known at the surface r = r^ and f.^^ is known at 
positions far away from the spacecraft. Since we are looking for equilibrium 
.solutions of Vlasov's equation, the distribution functions satisfy 



along particle trajectories. The object then is to calculate the trajectories of a 
selection of particles. Since we know the distribution function f far from the 
satellitCi automatically we know the distribution, since f is constant along a 
given trajectory. 

'Die Parker-Whipple® Inslde-oUt scheme makes use of this fact. Trajectories 
are initiated at the spacecraft surface and are traced backivards through the 
potential field to distant points where the distribution, f(v) is known. One advantage 
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6f this sehehne is that orbits which dc hot h&ve- ah endjjblht Oft thfe Spacecraft sUr- 
face are ai^otded. The neglect of such orbits is strictly >/alid only If the volume 
sp&c^bhdrge had no Ihflaonce on the pi^irtiOle tfdjefctories. 

It Should be emphaslaed that the machinery is contained In oUr numerical 
technique for calculating ambient Charge densities by COftstructtrtg distribution 
functions, f, in each spatial zone and taking its aeroth moment las opposed to j. 
Which is a first moment) 

P - Jt dV 

However, presently, we do not calculate this term, based upon the arguments 
presented in Section 2. 

3.2- -PotMitial Calculation 

In calculating the potential In three dimensions around an arbitrary object, 
a gridded method must be employed since the Specification of the surface is far 
tOo general for analytical or multipole techniques. Since satellites are the order 
of meters in length. We need at least 10 cm resolution as an upper bound in the 
Vicinity of the spacecraft. However, for determining particle orbits, the fields 
hundreds of meters away must also be known. In order to keep Storage down to 
a reasonable level, some type of variable griddlng must be employed. This 
precludes the use of any straightforward Fourier transform technique. One 
technique for achieving high resolution in the region around the object and still 
being able to handle vast quantities of space is through local mesh refinement. 
Finite difference approaches, however, have difficulty in meSh transition regions, 
especially when grid lines are terminated, and generally loSe an order Of accuracy 
in such regions. 

As a result of this, we decided upon a finite element approach using right 
parallelepiped elements and blended linear univariate edge interpolants. This 
permits the same degree of accuracy over the entire mesh, even though the mesh 
elements differ in size. It results in the standard trllinear Intei-polatlon Scheme 
for each element. 

The fundamental approach is. to solve Poisson'S equation 

by solving the associated variational principle 

j /dV l(W)^ + 4trp« I + /dS- (4»d+yd)^j (gj 
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The flrsWehto-lh the Ihteghahd cotrespdttds to the Laplhtlfiih ot>erat6r. The 
second term Id the volume spaoecHahge coftlrlbutioh. The Remaining terms are 
surface cohtributions. iJefehrlng to the surface charge and electric field, res“ 
pefctiVely. 

In the variational calculation, «e use locally defined basis sets, that Is, 
trlllnear InterpOlants within each cube-llke element. Since the finite element 
equations are derived from Eq. (2), different meSh volumes automatically receive 
the correct variational Weight. This ensures the maintenance of accuracy tlirough 
mesh transition regions. The problem of local mesh refinement IS approached by 
having grids within grids, that Is, a Chinese doll-like hierarchy of grids shown 
schematically In Figure 1. The theory of this technique is discussed In Blrkhoff 
etal,'^ and Cavendish. ® Imorder to have high computation speed, the linear 
equations resulting from the Variational principle (Eq. (2)) Ih the interface region 
were coded Up explicitly In a series of thirteen subroutines. These Same routines 
are used for interfacing any pair of the meshes. 



Figure 1. Cross-Section Of Grid, Showing ElrSt Four Embedded Meshes 
4. SAMPLE CALCULATIONS 

To demonstrate the capabilities 6f our 3D mOdel, We have performed tWo 
Sample calculations. First, we calculated the capacttance, Surface charge 
distribution and electric fields around a geometrically complex, conducting satel- 
llfcllke object. The object is shown in Figure 2, and the probiSiii Was gridded as 


326 



Figure 2. 3-D Modfel Spacecraft for Cap acitahce Ci»lMitn>ion 


shown in Fi^re L While it is olfectrically simple, being an equlpolential surface, 
it has sufflcieht geometrical complexity to demonstrate some of the features of 
our TRlLlN modeL The overall length of the object is 6 m. with 20 Cm resolutioh 
on the Surface. The outermost grid is 51 m long, and there are about 30. OOO 
variables in the problem. The outermost mesh had monopole 0 = S) boundary 
cohdltions imposed. Using an SOR routine, this problem took lessttiaiud min-to^ 
solve On the CDC 7600 at KirtlSnd Air Rorce Base. 

The capacitance calculated for this object is 03 pF. The surface area of the 
object is almost four times aS great as that of a Sphere of equivalent Capacitance 
(r = 75 cm). If placed ih an environment With a Charging current Of 10"9 A/cm^. 
this satellite-like object would charge to 10 keV in about 3 msec. The charge 
distribution is nonuhiform. as expected, with most of the charge on the panels 
Which have only 58 percent Of the Surface area. With the satellite charged to 
10 kV. the total charge oh the surface is approximately 2500 esu (-0. 9 nC). The 
average normal electric fields ort each panel in s'Uch a problem is —37. 5 V/cm 
While, on the body, it ranges from 2oV/cmto-40V/cm. 

The second sample calculation is of a simpler geometry, but has considerable 
physical complexity. Tlie object consists of a conducting cube. 60 cm on an edge, 
partially covered with a 1 mm insulation skin of dielectric constant unity. 

Figure 3 shows a picture Of the object. The object is placed in a 10 keV, 
hg » 10 cm-3 electron pUsmS with an assumed neutralising background. By the 
backward trajectory technique described above, incident electron currents on the 
object are determined. Charge impinging upon the dielectric skin is assumed to 
stick while charge landing Oh the exposed conducting surfaces is allowed to 
distribute itself in order to maintain the conductor aS ah fequtpotential surface, 
the potential on dielectric surfaces is related to that ort the conductor by the line 
integral Of the electric field through the surface. To add asymmetry and cause 
differential charging, we assume a light source at some large distance alohg the 
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poaitlve x«axl8* tliid light glvas rl88 to a 
photdcurriant of 4 nA/dm^ whicluactst4a disu- 
chai*ge that surface. 

initially the potential on the suj^face of 
the satellite rises very t^ulekly. Thia cor- 
reapdhds to surface charging on a timescale 
determined by the net capacitance to infinity. 
However, after about 50 msec,, the differen- 
tial charging of the conductor and the front 
and back dielectric surfaces dominate the 
calculation. T?he potential at three loca- 



tions as a function of time is plotted in 
Figure 4. We notice how the surface 
dielectric continues to charge, albeit at 


Figure 3. Spacecraft for Dynamicai 

Calculation, ^nly dashed area is 
bare metal, thv est Of the object is 
covered with a dielectric film 




Time (iee) 

Figiirc 4. Comparison Of Charge Buildup at liiree Different Sections of 
Illuminated Spacecraft 
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a v6ry slow rate, while the front surface dielectric discharges Substantially. 

T^he CbAductor also discharged, but more slowly than the illuminated dielectric. 
Figure 5 showS a potential contour map through the 3t, y plane. We can see that 
the Conductor is more than one thousand volts negative with respect to the front 
surface dielectric, while it is only a few hundred Volts positive With respect to 
the rear surface dark dielectric. This implies that the conductors Surface charge 
under the illuminated dielectric is of negative sign while underneath the rear 
dielectric the conductor's surface charge is of positive Sign. The potential dif- 
ference between front and rear dielectrics iS almost tWo kilovolts. Fields in the 
front dielectric are greater than 10^ Volts/cm. 



Figure 5. Potential Contour Plot Near the Space- 
craft after 0. 21 sec. Sunlight is incident from 
the right (jc- direction) 


5. CONCLUSIONS 

The prediction of surface potentials on complex satellites is a formidable 
task. Material properties, geometrical effects, ambient plasma, and photosheath 
space charge all play roles in determining surface potential distributions. 
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However, for the range of plasma parameters fl^eqUently found In Wagnetospherl. 
subStorms. vt lS justifiable to neglect the self- consistent ambient plasma space 
charge. This assumption permits the calculation of potentials In asymmetric 
three-dimensional geometries. The resultant calculations demonstrate such 
effects as net object charging with respect to infinity, differential chargi^. a 
charge redistribution on conductors. These first calculations presented here 
employ large Simplifications with reSpect to material properties. lOn ^ 

etc However, they Show that the concept Of three-dimensional Spacecraft charg- 
ing calculations iS a practical one. 
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